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ABSTRACT 

Context. The chemical composition of a molecular cloud changes dramatically as it collapses to form a low-mass protostar and cir- 
cumstellar disk. Two-dimensional (2D) chemodynamical models are required to properly study this process. 

Aims. The goal of this work is to follow, for the first time, the chemical evolution in two dimensions all the way from a pre-stellar 
core into a circumstellar disk. Of special interest is the question whether the chemical composition of the disk is a result of chemical 
processing during the collapse phase, or whether it is determined by in situ processing after the disk has formed. 
Methods. Our model combines a semi-analytical method to get 2D axisymmetric density and velocity structures with detailed ra- 
diative transfer calculations to get temperature profiles and UV fluxes. Material is followed in from the core to the disk and a full 
gas-phase chemistry network - including freeze-out onto and evaporation from cold dust grains - is evolved along these trajectories. 
The abundances thus obtained are compared to the results from a static disk model and to observations of comets. 
Results. The chemistry during the collapse phase is dominated by a few key processes, such as the evaporation of CO or the photodis- 
sociation of H 2 0. Depending on the physical conditions encountered along specific trajectories, some of these processes are absent. 
At the end of the collapse phase, the disk can be divided into zones with different chemical histories. The disk is not in chemical equi- 
librium at the end of the collapse, so care must be taken when choosing the initial abundances for stand-alone disk chemistry models. 
Our model results imply that comets must be formed from material with different chemical histories: some of it is strongly processed, 
some of it remains pristine. Variations between individual comets are possible if they formed at different positions or different times 
in the solar nebula. 

Key words, astrochemistry - stars: formation - circumstellar matter - planetary systems: protoplanetary disks - molecular processes 



"xt" 1 . Introduction was found to evaporate during the infall phase and to freeze out 

again in those parts of the disk that are colder than ~20 K. The 

' I , The formation of a low-mass protostar out of a cold molecu- higher binding energy of H2O keeps it in solid form at all times, 

■ lar cloud is accompanied by large-scale changes in the chemical except within ~10 AU of the star. Based on the time that the 

' composition of the constituent gas and dust. Pre-stellar cloud infalling material spends at dust temperatures between 20 and 

^ — I , cores are cold (~10 K), moderately dense (~10 4 -10 6 cm -3 ), 40 K, first-generation complex organic species were predicted to 

,— 1 ■ and irradiated onl y by the ambient i nterstellar radiation field fo rm abundantly on the gra in surfaces accord ing to the scenario 

^! ' (IShuetal.1 Il987t IpTFrancescoetall l2007t IBergin & Tafallal of lGarrod & Herbstl d2006l) and lGarrod et all d2008l) . 

• . |2007|). The material heats up as the core collapses, and the in- The current paper extends the chemical analysis to a full gas- 

P<j ; ner few hundred AU flatten out to form a circumstellar disk phase network, including freeze-out onto and evaporation from 

H , with densities that are orders o f magnitude higher than in the dust grains, as well as basic grain-surface hydrogenation reac- 

parent cloud dDullemond etalj |2007|). Meanwhile, the proto- t i ons . Combining semi-analytical density and velocity structures 

star infuses the disk and remnant envelope with large fluxes with detailed temperature profiles from full radiative transfer 

of ultraviolet and X-ray photons. The chemical changes aris- calculations, our aim is to bridge the gap between ID chemical 

ing from these evolving physical conditions ha ve been anal- mo dels of collapsing cores and 1 + ID or 2D chemic al models of 

ysed with one-dimensional spherical models (e.g.,|Aikawa etap T Tauri and Herbig Ae/Be disks dBergin et alJl2007l) . One of the 

|2008|). However, two-dimensional models are required to prop- k ey questions is whether the chemical composition of such disks 

erly describe the formation of the circumstellar disk and to anal- resu it s mainly from chemical processing during the collapse or 

yse the full chemical evolution from core to disk. whether it is determined by in situ processing after the disk has 

This paper is the second in a series of publications aiming formed, 

to model the chemical evolution from pre-stellar cores to cir- As reviewed by iDi Francesco et a D d2007h and 

cums tellar disks in two dimensions. The first paper dVisser et alJ IBergin &~T afalla (2007]), the chemistry of pre-stellar cores 

2009b, hereafter Paper I) contained a detailed description of our is well understood. Because of the low temperatures and 

semi-analytical model and an analysis of the gas and ice abun- the moderately high densities, many molecules are depleted 

dances of carbon monoxide (CO) and water (H2O). Most CO from the gas by freezing out onto the cold dust grains. The 
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main ice constituent is H2O, showing abundances of ~10 



ving 

relative to gas-phase H2 dTielens et alJI 199 it iPontoppidan et alJ 
20041). Other abund ant ices are C0 2 and CO dGibb euiil 2004 



1. Set initial conditions: 



Oberg et al] 1201 lal) . Correspondingly, the observed gas-phase 
abundances of H2O, CO and CO2 in pre - stellar cores are low 
(ISnell et all l2000l: iBacmann et all l200l ICaselli et all l20loh . 
Nitrogen-bearing species like N2 and NH3 are g enerally less 
depleted than carbon- and oxygen-bearing ones (Tafall a et alJ 
2004), partly because they require a longer time to be formed 
in t he gas and therefore have not yet had a chan ce to freeze 
out dAikawa et alJ 120011: iDi Francesco et alJ 120071) . and partly 
because CO is no longer present as a destructive agent 
(iRawlings et alJfl992l) . The observed d epletion factors are well 
reproduced wit h ID chemical models (Bergin & Langerlll997t 
iLee et al . 200|). 

The collapse phase is initially characterised by a grad- 
ual warm-up of the material, resulting in the evaporation 
of the ices according to their respective bindi ng energies 



(Ivan Dishoeck & Blake! 1998Hj0rgensen et alJ2004h . The higher 
temperatures also drive a rich chemistry, especially if it gets 
warm enough to evapo rate H2O and orga nic species like 
CH3OH and HCOOCH3 (ICharnlev et al.lll992l) . Spherical mod- 
els of the chemical evolution during the co l lapse phase (e.g.. 



Rod gers & Charnleyl [20031: ILee et all 
2006t lAikawa et alJl2008t iGarrod et al 



2004; IGarrod & H erbst 



2008) are successful at 



explaining the observed abundances at scales of several thou- 
sand AU, where the envelope is still close to spherically sym- 
metric, but they cannot make the transition from the ID spheri- 
cally symmetri c envelope to the 2D axi symmetric circumstellar 
disk. Recently. Ivan Weeren et al. I J2009¥ followed the chemical 
evolution within the framework of a 2D hydrodynamical sim- 
ulation and obtained a reasonable match with observations of 
protostars. Even though their primary focus was still on the en- 
velope and not on the disk, they did show how important it is to 
treat the chemical evolution during low-mass star formation in 
more than one dimension. 

Once the phase of active accretion from the envelope comes 
to an end, the circumstellar disk settles into the compara- 
tively static T Tauri or Herbig Ae/Be phase. So f ar only some 
simple molecules have b e en detected in disks (Dutrey et al. 
1997t iKastner etalJ 1 1991 iThi et al J l2004t IQberg et al] 120101 
201 lbl) . but the inventory is expected to grow now that the 
Atacama Large Millimeter Array (ALMA) is becoming oper- 
ational. Chemically, disks can be divided into three layers: a 
cold zone near the midplane, a warm molecular layer at inter- 
mediate altitudes, a nd a photon-dominated region at the surface 
(Berg in et al.l[20 07*). However, it is unknown to what extent the 
chemical composition of the disk is affected by the chemical evo- 
lution in the collapsing core at earlier times, or how much chem- 
ical processing has been experienced by material ending up in 
planetary and cometary building blocks. 

This paper aims to provide an important step towards filling 
in this missing link by following the chemical evolution all the 
way from a pre-stellar cloud core to a circumstellar disk in two 
spatial dimensions. The physical and chemical models are de- 
scribed in Sects.[2]and[3] The chemistry during the pre-collapse 
phase is briefly reviewed in Sect. [4] followed by an extensive 
discussion of the collapse-phase chemistry in Sect. [5] The abun- 
dances resulting from the collapse phase are compared to in situ 
processing in a static disk in Sect. [6] Sections|7]and[8]discuss the 
implications of our results for the origin of comets and address 
some caveats in the model. Conclusions are drawn in Sect. [9] 



2. Compute density and 
velocity structure: 
n(R,z,f), v(R,z,t) 




3. Compute temperature 
and radiation field: 
r d (R,z,f), F uv (R,z,f) 




Semi -analytical L----^*! 


Numerical simulation 


4. Compute infall 




5. Run chemistry for 


trajectories 




each parcel 



6. Transform parcel 
abundances to 2D profiles 



Fig. 1. Step-wise summary of our 2D axisymmetric collapse 
model. Steps 2 and 4 are semi-analytical, while steps 3 and 5 
consist of detailed numerical simulations. 



2. Collapse model 

2.1. Step-wise summary 

Our sem i-analytical collapse model is described in detail in 
Paper I and summarised in Fig.Q] The model starts with a singu- 
lar isothermal sphere characterised by a total mass Mo, an effec- 
tive sound speed c s , and a uniform rotation rate Qo. As soon as 
the collapse starts, at t = 0, the rotation causes the infalling mate- 
rial to be deflected towards the equatorial midplane. This breaks 
the spherical symmetry, so the entire model is run as a two- 
dimensional axisymmetric sys tem. The 2D density and velocity 
profile s fol low the solutions of [ Shul dl977l) . ICassen & Moosman 
(Il98lh and lTerebeveTall (1 19841 for an inside-out collapse with 
rotation. After the disk is first formed at the midplane, it evolves 
by ongoing accretion from the collapsing core and by viscous 
spreading to conserve angular mo mentum ( Shakura & Sunvaevl 
fl97llLvnden-Bell & Pringi3[T97l . 

Taking the 2D density profiles from step 2, an d adopting the 
appropriate size and lumino sity for the protostar dAdams & Shul 
1986: [Young & Evansl2005l) . the next step consists of computing 
the dust temperature at a number of time steps. This is done with 
the r adiative transfer code RADMC dDullemond & Dom inik 
12004 . which takes a 2D axisymmetric density profile but fol- 
lows photons in all three dimensions. RADMC also computes 
the full radiation spectrum at each point in the axisymmetric 
disk and remnant envelope, as required for the photon-driven re- 
actions in our chemical network (Sects. 12.31 and 13. U . The gas 
temperature is set equal to the dust temperature throughout the 
disk and the envelope. This is a poor assump tion in the surface of 
the disk and the inner parts of the envelope (Kam p & Dullemondl 
l2004t [fonkheid et alJl200l . the consequ ences of which are ad- 
dressed in Sect. [8] As argued in iPaper it the disk-envelope ac- 
cretion shock is not important for the chemical evolution of the 
material in the disk at the end of the collapse phase. It is there- 
fore not incorporated in our model. 

Given the dynamical nature of the collap se, it is easiest to 
solve the chemistry in a Lagrangian frame. In P aper 1 the enve- 
lope was populated with several thousand parcels at t = and 
these were followed in towards the disk and star. The method is 
reversed here: a regular grid of parcels is defined at the end of 
the collapse and the parcels are followed backwards in time to 
their position at t = 0. The reason for doing so is that the 2D 
abundance profiles at the end of the collapse, when the disk is 
fully formed, are of greater interest than those at the onset of 
collapse. In either case, step 4 of the model produces a set of 
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infall trajectories with densities, temperatures and UV intensi- 
ties as a function of time and position. These data are required 
for the next step: solving the time-dependent chemistry for each 
individual parcel. Although the parcels are followed backwards 
in time to get their trajectories, the chemistry is computed in the 
normal forward direction. The last step of our model consists of 
transforming the abundances from the individual parcels back 
into 2D axisymmetric profiles at any time step of interest. 

In lPaper j , the model was run for a grid of initial conditions. 
In the current paper, the analysis is limited to our standard set 
of parameters: M = 1.0 M , c s = 0.26 km s -1 and Qp = 10~ 14 
s" 1 . T he viscosity parameter a is fixed at 0.01 (lHartmann et al.l 
1998). Section □ contains a brief discussion on how the results 
may change for other parameter values. 

2.2. Differences with Paper J 

The current version of the model contains several improve- 
ments over the version use d in IPaper iL described in detail by 
IVisser & Dullemondl d2010b . Most importantly, the model now 
correctly treats the problem of sub-Keplerian accretion onto a 
2D disk. Material falling onto the disk along an elliptic orbit has 
sub-Keplerian angular momentum, so it exerts a torque on the 
disk that r esults in an inward push. Several solutions are avail - 
able (e.g., ICassen & Moosmanl fT98lt iHueso & Guilloll l2005l) . 
but these are not suitable for our 2D model. The ad-hoc solu- 
tion from Paper I provided the appropriate qualitative physical 
correction - increasing the inward radial velocity of the disk ma- 
terial - but it did not properly conserve angular momentum. The 
updated model uses a new, fully consistent solution, derived di- 
rectly from the equ ations for the conservat ion of mass and an- 
gular momentum d Visser & Dullemondl2010l) . It results in disks 
that are typically a factor of a few smaller than those obtained 
previously, but the disk masses are mostly unchanged. The new 
disks are a few degrees colder in the inner part and warmer in 
the outer part, which may further affect the chemistry. 

Other changes to the model include the definition of the 
disk-en velope boundary and the shape of the outflow cavity. In 
Paper I, the disk-envelope boundary was defined as the surface 
where the density of the infaUing envelope material equals that 
of the disk. The model now uses the surface where the ram pres- 
sure of the infaUing material equ als the thermal pressure of the 
disk (Visser & Dullemond 2010), providing a more physically 
correct description of w here material becom es part of the disk. 
The collapse solution of Terebev et al.l (1 19841) does not produce 
an outflow cavity, so one is added manually. The original model 
used a cavity with straight walls, but observations and theoret- 
ical predictions show that curved walls are more appropriate 
(IVelusamv & Langerlfl998l: ICanto et al.ll2008l) . The adopted de- 
scription of the outflow wall is now 

z=(0 - 191AU) (i)"W' s ■ m 

with R and z in cylindrical coordinates and f acc = Mq /M the time 
required for the entire envelope to accrete onto the star and disk. 
The f -3 dependence is chosen so that the outflow starts narrow 
and becomes wider as the collapse proceeds. The full opening 
angle at r acc is 33.6° at z = 1000 AU and 15.9° at z = 10 000 AU. 
There are some streamlines in the infaUing envelope that hit the 
cavity wall, but the mass flowing along these lines is insignificant 
(~1%) compared to the mass accreting onto the star and the disk. 



2.3. Radiation field 

Photodissociation and photoionisation by ultraviolet (UV) radi- 
ation are important processes in the hot inner core and in the 
surface layers of the disk . The tem perature and luminosity of the 
protostar change in time (Paper I), so neither the strength nor the 
spectral shape of the radiation field is constant. In addition, the 
spectral shape changes as the radiation passes through the disk 
and remnant envelope. 

The most accurate way of obtaining the time- and location- 
dependent photorates is to multiply the cross section for each re- 
action by the UV field at each grid point. The latter can be com- 
puted from 2D radiative transfer at high spectral resolution. As 
this is too computationally demanding, several approximations 
have to be made. The first one is to assume that the wavelength- 
dependent attenuation of the radiation field by the dust in the 
disk and envel ope can be represented by a single factor y for 
each reaction (van Dishoe ck et a 1. 2006). The rate coefficient for 
a given photoreaction at spatial coordinates r and 9 and at time t 
can then be expressed as 

V' e > f ) = fe Ph(£) ex P(-r^v), (2) 

with Ay the visual extinction (see below). The unshielded rate 
coefficient is calculated at the stellar surface (&* h ) by multiplying 
the cross section of the reaction by the blackbody flux at the ef- 
fective temperature of the protostar, which va ries betw een 4500 
and 5500 K for the duration of the collapse (IPaper J) . Any ex- 
cess UV from the disk-star boundary is not included. The term 
{r/R„y 2 accounts for the geometrical dilution of the radiation 
from the star across a distance r. The factor y is discussed in 
Sect. ED 

In order to apply Eq. (0, the extinction towards each point 
(r, 0) is calculated with RADMC at a low spectral resolution of 
one frequency point per eV. The local UV flux from RADMC at 
(r, 0), Fradmc (integrated from 6 to 13.6 eV), can be considered 
the product of the unattenuated flux and an extinction term, 

^RADMcCn 0, t) = F exp(-T UV , e ff) = F„ j exp(-T UV , e ff) ■ 

(3) 

F* is the time-dependent UV flux at the stellar surface. The ef- 
fective UV optical depth, Tuv.eff, is an average over the many 
possible paths a photon can follow from the source to the point 
(r, ff). The large number of photons propagated through the grid 
by RADMC (typically 10 5 ) ensures a statistical sampling of all 
possible trajectories. The UV optical depth is converted to the 
visual extincti on Ay through t he standard relationship Ay = 
ruv.eff/3.02 dBohlin et al.lll978l) . 

3. Chemical network 

The basis of our ch emical network is the UMIST06 dat abase 
(IWoodall et alJl2007l) as modified bv lBradereretail (120091) . ex- 
cept that X-ray chemistry is not included. The modified net- 
work includ es charge-exchange re actions between positive ions 
and grains (Malon ev et alJ fl996l) . which can be important in 
the disk's inner few AU. T he cosmic-ray i onisation rate of 
H 2 is set to 5 x 10~ 17 s _1 (Dalgarnol 120061) . Attenuation of 
cosmic rays at column densi ties of more than ~100 g cm -2 
(lUmebavashi & Nakandll98ll) is not taken into account, even 
though this can become relevant for the midplane chemistry 
in the inner few AU of massive disks. The network contains 
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Table 1. Elemental abundances: x(X) = n(X)/nn. 



Species 


Abundance 


Species 


Abundance 


H 2 


5.00C-1) 


Si 


2.74(-9) 


He 


9.75(-2) 


P 


2.16(-10) 


C 


7.86(-5) 


S 


9.14(-8) 


N 


2.47(-5) 


CI 


1.00(-9) 





1.80(-4) 


At 


3.80(-8) 


Nc 


1.40(-6) 


Ca 


2.20(-8) 


Na 


2.25(-9) 


Cr 


4.90(-9) 


Mg 


1.09(-8) 


Fc 


2.74(-9) 


Al 


3.10(-8) 


Ni 


1.80(-8) 



Notes. The notation a(b) means a x lO 6 



162 neutral species, 251 cations and six anions, built up out 
of 18 elements. The initial composition is fully atomic, except 
that hydrogen start s as Hg . Elemental abundances are adopted 
from Aikawa etal. ( 20081) . Additional values are taken from 
iBruderer et al. " 3200S ), but reduced by a factor of 100 to have 
a consistent set of low-metal abundances. Table Q] lists the ele- 
mental abundances relative to the total hydrogen nucleus den- 
sity: n H = «(H) + 2n(H 2 ). 

In order to set the chemical composition at the onset of col- 
lapse (f = 0), the initially atomic gas is evolved for a period of 1 



Myrat«H = 8xl0 4 cirT' > and Ft, = = 10 K. The visual extinc- 
tion during this pre-collapse phase is set to 100 mag to disable 
all photoprocesses, except for a minor contribution from cosmic- 
ray-induced photons. The resulting solid and gas-phase abun- 
danc es are consistent with th ose observed in pre-stellar cores 
fe.g jDi Francesco et al.ll2007h . They form the initial conditions 
for the collapse phase for all infalling parcels. In the remainder 
of this paper, t — always refers to the onset of collapse, follow- 
ing the 1 Myr pre-collapse phase here described. 

3.1. Photodissociation and photoionisation 

Photodissociation and photoionisation by UV radiation are im- 
portant processes in the inner disk and inner envelope. Their 
rates are given by Eq. (0, using the effective extinction from 
Eq. (0. The unshielded rates at the stellar surface (k* h ) are cal- 
culated for a blackbody spectrum at the time-dependent stel- 
lar temperatures from iPaper it The required cross sections are 
taken from our freely available database^] V alues for the extinc- 
tion fa ctor y from Eq. (0 were tabulated by van Dishoec k et alj 
(2006) for three different types of radiation fields, of which the 
4000 K blackbody is most appropriate for our model (Sect. l2~3l >. 
The final output of this procedure is a 3D array (two spatial 
coordinates and time) with rate coefficients (£ p h) for each pho- 
toreaction. When computing the infall trajectories for individual 
parcels (step 4 in Fig.0, the rate coefficients at all points along 
each trajectory are obtained from a linear interpolation. 

The photodissociation of H2 and CO requires some special 
treatment. Both processes occur exclusively through discrete ab- 
sorption lines, so self-shielding plays an important role. The 
amount of shielding for H2 is a function of the H2 column den- 
sity; for CO, it is a function of both the CO and the H2 column 
density, because some CO lines are shielded by H2 lines. The 
effective UV extinction from Eq. (0 can be converted into a 
total hydrogen column de nsity through Ay = Tijy. e ff /3.02 and 
N H = 1.59 x 10 21 A V cm- 2 dBohlin et alJI 1 978t iDiplas & Savagd 
119941) . Most hydrogen along each photon path is in molecu- 



Table 2. Pre-exponential factors and binding energies for se- 
lected species in our network. 



Species 


v(s-') 


Eb/* (K) 


Reference 


H 2 


1 x 10 15 


5773 


Fraseretal. (2001) 


CO 


9 x 10" 


855 


Bisschop et al. (2006) 


N 2 


1 x 10" 


800 


Bisschop et al. ('2006) 


O, 


9x 10" 


912 


Acharvva et al. (2007) 



lar form, so the effective N(H2) towards each spati a l grid point 
is ~0.5/V H - Equation (37) from iDraine & Bertoldil U996) then 
gives the amount of self-shielding for H2. The unshielded disso- 
ciatio n rate is computed a ccording to the one-line approximation 
from van Dish oeck] (1 19871). scaled so that the rate is 4.5 x 10" 
s _1 in the standard Draine (1978) field. The photodissociation 
rate of CO is c omputed from our new shielding functions and 
cross sections (I Visser et al.ll2009ah . Effective CO column den- 
sities for the self-shielding calculation are derived from the H2 
column densities by assuming an average N(CO)/N(H2) ratio of 
10~ 5 , accounting for the partial dissociation of CO along each 
photon path. The results are not sensitive to the e xact choice of 
average CO abundance (van Zadelhoff etal . 200i. 



3.2. Gas-grain interactions 

All neutral species other than H, H2, He, Ne and Ar are allowe d 
to freeze out onto the dust according to ICharnlev et afl (120011) . 
In cold, dense environments - such as our model cores before 
the onset of collapse - obser vations show H2O , CO and CO2 to 
be the most abundant ices (Boogert et ai] l201ll) . As the temper- 
ature rises during the collapse, the ices evaporate according to 
their respective binding energies. However, the presence of non- 
volatile species like H2O prevents the more volatile species like 
CO and CO2 from evap orating entirely: some CO and CO? gets 
trapp ed in the H2O ice dHasegawa & Herbsllll993tlFav olle et al. 
1201 ll) . The results from Paper!] suggested that this is required to 
explain the presence of CO in solar-system comets. 

Incorpora ting ice trapping in a chemical ne twork is a non- 
trivial task dViti et alJl2004t iFavolle et a"Dl201ll) . so the process 
is ignored for now. Instead, desorption of all species is treated 
according to the zeroth-order rate equation 



#thdes(X) = 4nain f(X)v(X)N ss cxp 



E b (X) 



kT d 



(4) 



where Ta is the dust temperature, a gl = 0.1 /mi the typical 
grain radius, and n gr = 1 x 10~ 12 «h the grain number density. 
The canonical pre-e xponential factor, v, for first-ord er desorp- 
tion is 2 x 10 12 s" 1 dSandford & Allamandolalll993h . which is 
used for all ices except the four listed in Table [2] In order to get 
the zeroth-order pre-exponential factor, v is multiplied by N ss , 
the number of binding sites per unit grain surface. Assuming 



JV SS 



8 x 10 14 cm z , the number of binding sites per grain 



is Nb = 1 x 10 6 for our 0.1 yum grains. The binding energies 
of species o ther than the four in Tabl e [2] are set to the values 
tabulated by Sandford & Allamandola (1993) and Aikawa et al. 
(1 19971) . Species for which the binding energy is unknown are as- 
signed the binding energy and the pre-exponential factor of H2O. 
The dimensionless factor / in Eq. (01 ensures that each species 
desorbs according to its solid-phase abundance, and changes the 
overall desorption behaviour from zeroth to first order when 
there is less than one monolayer of ice: 



http://www.strw.leidenuniv.nl/~ewine/photo 



f(X) = 



n s (X) 



max(n ice ,/VbHgr) 



(5) 
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with « S (X) the number density of ice species X, Nb = 10 6 the typ- 
ical number of binding sites per grain, and «j ce the total number 
density (per unit volume of cloud or disk) of all ice species com- 
bined. Section [8]briefiy discusses how our results would change 
if trapping were included. 

In addition to thermal desorption, our model includes des- 
orption induced by UV photons. Laboratory experiments on the 
photodesorption of H2O, CO and CO2 all produce a yield of Y sa 
10" 3 molecu les per grain per incident UV photon (Ober get alj 
120071 l2009al). while the y ield for N 2 is an order of magnitude 
lower (Oberget al. 2009b). For all other ice species in our net- 
work, whose photodesorption yields have not yet been deter- 
mined experimentally or theoretically, the yield is set to 10" 3 . 
The photodesorption rate is 



AphdesPO = na„n f(X)Y(X)F exp(-r uv ,eff) , 



(6) 



with / the same factor as for thermal desorption. The unattenu- 
ated UV flux (F^) and the effective UV extinction (ruv.eff) follow 
from Eq. (|3). Photodesorption occurs even in strongly shielded 
regions because of cosmic-ray-induced photons, which is ac- 
counted for by settin g a lower limit of 10 4 cirT 2 s _1 to the prod- 
uct F exp(-Tuv,eff ) CSh en et al.ll2004 . 

The chemical reactions in our model are not limited to the 
gas phase. As us ual, the network includes the grain-surface 
formation of H 2 dBlack & van Dishoeckl 1 987[). I nspired by 
Bergin & Langerl (fl997) andlHollenbac h et al. (2009), it also in- 
cludes the grain-surface hydrogenation of C to CH4, N to NH3, 
O to H2O, and S to H2S. The hydrogenation is done one H atom 
at a time and is always in competition with thermal and photon- 
induced desorption. The formation of CH4, NH3, H2O and H2S 
does not have to start with the respective atom freezing out. For 
instance, CH freezing out from the gas is also subject to hydro- 
genation on the grain surface. The rate of each hydrogenation 
step is taken to be the adsorption rate of H from the gas multi- 
plied by the probability that the H atom encounters the atom or 
molecule X to hydrogenate: 



8kT B 



(7) 



with To the gas temperature. The factor /' serves a similar pur- 
pose as the factor / in Eqs. © and ((6]). Since the hydrogenation 
is assumed to be near-instantaneous as soon as the H atom meets 
X before X desorbs, X is assumed to reside always near the top 
layer of the ice. Hence, the abundance of solid X should not be 
compared to the total amount of ice (as in /), but to the combined 
abundance of other ice species that can be hydrogenated: 



f'(X) = 



n s (X) 



max(n hy di-o, iVb«gr) 



(8) 



with «hydro the sum of the solid abundances of the eleven species 
subject to hydrogenation: C, CH, CH 2 , CH 3 , N, NH, NH 2 , O, 
OH, S and SH. The main effect of this hydrogenation scheme is 
to build up an ice mixture of simple saturated molecules dur- 
ing the pre-collapse phase, as is f ound observationally (e.g., 
Tiele ns et alJll99lHGibb etafll2004l) . 

In reality, grain-surface chemistry is not limited to the sim- 
ple hydrogenation steps included here. For example, CO can 
react with H to form H 2 CO and CH3OH or with OH to form 
C0 2 dFuchs et alJ l2009[ IIoppoIo et al l 1201 ll) . and 2 can be 
hydrogenated to H2O (lloppolo et al.l 120081) . Grains also play 
an important role in the formation of mo re complex species 
dGarrod & Herbstll2006t iGarrod et alJ l2008). However, none of 



these reactions can be implemented as easily as the hydrogena- 
tion of C, N, O and S. In addition, the main focus of this pa- 
per is on simple molecules whose abundances can be well ex- 
plained with conventional gas-phase chemistry. The effects of 
more complex grain-surface chemistry will be explored in detail 
in a follow-up paper. 



4. Results from the pre-collapse phase 

This section, together with the next two, contains the results 
from the gas-phase chemistry in our collapse model. The chem- 
istry during the pre-collapse phase is briefly discussed first. The 
chemistry during the collapse is analysed in detail for one par- 
ticular parcel in Sect. 15. II and then generalised to others in Sect. 
15.21 Finally, the collapse chemistry is compared to a static disk 
model in Sect. [6] The results in the current section are all consis- 
tent with available observatio nal constraints on pre-stellar cores 
(e.g. JDi Francesco et al.ll2007h . 

During the 1.0 Myr pre-collapse phase, most of the initial 
atomic oxygen freezes out and is hydrogenated to H2O ice. 
Meanwhile, H2 is ionised by cosmic rays. The resulting Hi re- 
acts with H2 to give H3 , which reacts with the remaining atomic 
O to ultimately form 2 : 



O + H+ 
OH + + H 2 
H 2 + + H 2 
H 3 + + e 

OH + O 



OH + + H 2 , 
H 2 + + H, 
H 3 + + H, 
OH + H 2 /2H, 
2 + H. 



(9) 
(10) 

(ID 
(12) 
(13) 



The O2 thus produced freezes out for the most part. At the onset 
of collapse, the four major oxygen reservoirs are H2O ice (44%), 
CO ice (34%), 2 ice (16%) and NO ice (3%). 

The oxygen chemistry is tied closely to the carbon chemistry 
through CO. It is initially formed in the gas phase from CH2, 
which in turn is formed from atomic C: 



C + H 2 
CH 2 + O 



CH 2 + hv, 
CO + 2H. 



(14) 
(15) 



Another early pathway from C to CO is powered by HJ and goes 
through an HCO + intermediate: 



C + Hi 



CH + 
CH} 



H 2 
H, 



CH + 
CH 2 + 
CHi 



H 2 , 

H, 

H, 



ch; + o 

HCO + + C 



HCO + + H 2 , 
CO + CH + . 



(16) 
(17) 
(18) 
(19) 
(20) 



The formation of CO through these two pathways accounts for 
most of the pre-collapse processing of carbon: at t = 0, 82% of 
all carbon has been converted into CO, of which 97% has frozen 
out onto the grains. Most of the remaining carbon is present as 
CH4 ice (14% of all C), formed from the rapid grain-surface 
hydrogenation of atomic C. 

The initial nitrogen chemistry consists mostly of converting 
atomic N into NH3, N2 and NO. The first of these is formed 
on the grains after freeze-out of N, in the same way that H2O 
and CH4 are formed from adsorbed O and C. The model re- 
sults show two pathways leading to N 2 . The first one starts 
with the cosmic -ray dissociation of H2 into H + and H", which 
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ICravens & Dalgarnd(ll978l) estimated to have a 0.015% chance 
of happening for each cosmic -ray ionisation event: 



H 2 + f 
N + H" 
NH + N 



H + + H" 
NH + e 
N 2 + H. 



(21) 
(22) 
(23) 



The other pathway couples the nitrogen chemistry to the carbon 
chemistry. It starts with Reactions (TT~6b— (fl~8b to form CH3 , fol- 
lowed by 



CHI + e 



CH 



H 2 /2H, 
CH + N — > CN + H, 
CN + N -> N 2 + C. 



(24) 
(25) 
(26) 



The nitrogen chemistry is also tied to the oxygen chemistry, 
forming NO out of N and OH: 



N + OH -> NO + H, 



(27) 



with OH formed by Reaction (T2]i. Nearly all of the N 2 and NO 
formed during the pre-collapse phase freezes out. At t — 0, solid 
N 2 , solid NH 3 and solid NO account for 41, 32 and 22% of all 
nitrogen. 

5. Results from the collapse phase 

The collapse-pha se chem istry is run for the standard set of model 



parameters from Paper I: Mo = 1 .0 M , c s 
Qo = 10~ 14 s _1 . The core collapses in f acc 



0.26 km s and 
2.52 x 10 5 yr to 



form a 0.76 M star surrounded by a compact disk (47 AU outer 
radius) of 0.13 M . At t — f acc , the ra dius at which H 2 Q freezes 
out (the snowline) lies at 5.1 AU. The iToomrd (fl 964) Q param- 
eter ranges from 130 at 0.1 AU to a minimum of 1.4 at 28 AU, 
so the disk remains gravitationally stable by a small margin. 

5. 1 . One single parcel 

The chemical evolution is first discussed in detail for one par- 
ticular infalling parcel of material. It starts near the edge of the 
cloud core, at r — 6710 AU from the center and 8 = 48.8° de- 
grees from the z axis (R = 5050 AU, z = 4420 AU). Its trajectory 
terminates at t — f acc in the inner part of the disk, at R = 6.3 AU 
and z = 2.4 AU, about 0.2 AU below the surface. The phys- 
ical conditions encountered along the trajectory are plotted in 
Fig. |2 The radiation field is characterised by the unattenuated 
flux and the visual extinction, where the former is expressed 
as a scaling factor relative to the average flux in the interstel- 
lar mediu m (ISM): y = Fq/Fi$m, with Fism = 8 x 10 7 photons 



cm" 2 s 1 (lDraind f~1978). This figure also shows the abundances 



of the main oxygen-, carbon- and nitrogen-bearing species. The 
four panels in the right column are regular plots as function of 
R (the coordinate along the midplane), tracing the parcel's tra- 
jectory from t = (R = 5050 AU) to t = f acc (R = 6.3 AU). 
The infall velocity of the parcel increases as it gets closer to the 
star, so the physical conditions and chemical abundances change 
more rapidly at later times. Hence, the panels in the left column 
are plotted as a function of f acc - 1 : the time remaining before the 
end of the collapse phase. In each individual panel, the parcel 
essentially moves from right to left. 

Figure[3]presents a schematic overview of the parcel's chem- 
ical evolution. It shows the infall trajectory of the parcel and the 
abundances of several species at four points along the trajectory. 
The physical conditions and the key reactions controlling those 



abundances are also listed. Most abundance changes for individ- 
ual species are related to one specific chemical event, such as the 
evaporation of CO or the photodissociation of H 2 0. The remain- 
der of this subsection discusses the abundance profiles from Fig. 
[2] and explains them in the context of Fig. [3] 

5.1.1. Oxygen chemistry 

At the onset of collapse (f = 0), the main oxygen reservoir is 
solid H 2 at an abundance of 8 x 10~ 5 relative to «h (Sect. |4]i. 
The abundance remains constant until the parcel gets to point 
C in Fig. [3j where the temperature is high enough for the H 2 
ice to evaporate. The parcel is now located close to the outflow 
wall, so the stellar UV field is only weakly attenuated (Ay = 0.9 
mag). Hence, the evaporating H 2 is immediately photodissoci- 
ated into H and OH, which in turn is dissociated into O and a 
second H atom. At R = 17 AU (23 AU inside of point C), the 
dust temperature is 150 K and all solid H 2 is gone. Moving 
in further, the parcel enters the surface layers of the disk and is 
quickly shielded from the stellar radiation (Ay = 10-20 mag). 
The temperature decreases at the same time to 1 14 K, allowing 
some H 2 ice to reform. The final H 2 ice abundance at f acc 
(point D) is4x 10~ 8 . 

The dissociative recombination of HaO + (formed by 
Reactions (l9l-(fTTl>) initially maintains H 2 in the gas at an 
abundance of 7 x 10~ 8 at t = 0. Following the sharp increase 
in the overall gas density at t — 2.1 x 10 5 yr (Fig. [2j, the freeze- 
out rate increases and the gas-phase H 2 abundance goes down 
to 3 x 10~ 10 at point A in Fig. [3] Moving on towards point B, the 
evaporation of 2 from the grains enables a new H 2 formation 
route: 



2 + C + — > CO + o + , 
+ + H 2 -> OH + + H, 



(28) 
(29) 



followed by Reactions (TTQb and (fTTT i to give HsO + , which re- 
combines with an electron to give H 2 0. The H 2 abundance 
thus increases to 3 x 10~ 9 at R = 300 AU. Further in, at point 
C, solid H 2 comes off the grains as described above. However, 
photodissociation keeps the gas-phase abundance from growing 
higher than ~10~ 7 . Once all H 2 ice is gone at R = 17 AU, the 
gas-phase abundance can no longer be sustained at 10~ 7 and it 
drops to 3 x 10~ 12 . Some H 2 is eventually reformed as the par- 
cel gets into the disk and is shielded from the stellar radiation, 
producing a final abundance of 9 x 10~ 9 relative to ria- 

Another main oxygen reservoir at t = is 2 . Its solid and 
gas-phase abundances in the model are 1 x 10~ 5 and4x!0~ 7 , con- 
sistent with available observat ional constraints ( Goldsmit h et al.l 
1 1 985L l20TTt iFuente et al.l 1993h . Gaseous 2 gradually continues 
to freeze out until it reaches a minimum gas-phase abundance of 
3 x 10~ I() just inside of point A. The temperature at that time is 
1 8 K, enough for 2 to slowly start evaporating thermally. The 
gas-phase abundance is up by a factor of ten by the time the 
dust temperature reaches 23 K, about halfway between points A 
and B. The evaporation is 99% complete as the parcel reaches 
R = 460 AU, about 140 AU inside of point B. The gas-phase 
abundance remains stable at 1 x 10~ 5 for the next few hundred 
years. Then, as the parcel gets closer to the outflow wall and into 
a region of lower extinction, the photodissociation of 2 sets in 
and its abundance decreases to 2 x 10~ 8 at point C. The evapo- 
ration and photodissociation of H 2 at that point enhances the 
abundances of OH and O, which react with each other to replen- 
ish some 2 . As soon as all the H 2 ice is gone, this 2 produc- 
tion channel quickly disappears and the 2 abundance drops to 
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Fig. 2. Physical conditions («h, T^, x and Ay) and abundances of the main oxygen-, carbon- and nitrogen-bearing species for the 
single parcel from Sect. 15.11 as function of time before the end of the collapse (left) and as function of horizontal position during 
the collapse (right). The grey bars correspond to the points A, B and C from Fig. [3] In each panel, the parcel essentially moves from 
right to left. 
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Fig. 3. Overview of the chemistry along the infall trajectory (solid black curve) of the single parcel from Sect. 15.11 The solid and 
dashed grey lines denote the surface of the disk and the outflow wall, both at t — f acc = 2.52 x 10 5 yr. Physical conditions, abundances 
(black bars: gas; grey bars: ice) and key reactions are indicated at four points (A, B, C and D) along the trajectory. The key processes 
governing the overall chemistry at each point are listed in the bottom right. The type of each reaction is indicated by colour, as listed 
in the top left. 



1 x 10 . Finally, when the parcel enters the top of the disk, O2 
is no longer photodissociated and its abundance goes back up to 
4 x 10~ 8 at point D. 

The abundance of gas-phase OH starts at 3 x 10" 7 . Its main 
formation pathway is initially the dissociative recombination of 
H30 + (Reaction (fT2l), and its main destructors are O, N and Hi. 
The increase in total density at 2.1 x 10 5 yr speeds up the de- 
struction reactions, and the OH abundance drops to 3 x 10" 10 
at point B. The evaporation of solid OH then briefly increases 



the gas-phase abundance to 1 x 10 . When all of the OH has 
evaporated at R = 300 AU, the gas-phase abundance goes down 
again to 5 x 10~ 10 over the next 150 AU. As the parcel contin- 
ues towards and past point C, the OH abundance is boosted to 
a maximum of 1 x 10~ 6 by the photodissociation of H2O. The 
high abundance lasts only briefly, however. As the last of the 
H2O evaporates and gets photodissociated, OH can no longer be 
formed as efficiently, and it is dissociated itself. At the end of the 
collapse, the OH abundance is ~10~ 14 . 
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The fifth main oxygen-bearing species is atomic O itself. Its 
abundance is 7 x 1CT 7 at t — and 1 x 1CT 4 at t — ? acc , ac- 
counting for respectively 0.4 and 56% of the total amount of 
non-refractory oxygen. Starting from t — 0, the O abundance re- 
mains constant during the first 2.0 x 10 5 yr of the collapse phase. 
The increasing overall density then speeds up the reactions with 
OH (forming O2) and H3 (forming OH + ), as well as the adsorp- 
tion onto the grains, and the O abundance decreases to a low of 
2 x 10~ 8 just before point A. The abundance goes back up due 
to the evaporation of CO and, at point B, of O2 and NO, with O 



formed from the following reactions: 

CO + He + -> C + + O + He, (30) 

NO + N -> N 2 + O, (31) 

2 + CN -> OCN + O. (32) 



Heading on towards point C, the photodissociation of O2, NO 
and H2O further drives up the amount of atomic O to the afore- 
mentioned final abundance of 1 x 10~ 4 . 

5.1.2. Carbon chemistry 

With solid and gas-phase abundances of 6 x 10~ 5 and 2 x 10 -6 
relative to hh, CO is the main form of non-refractory carbon 
at the onset of collapse. CO is a very stable molecule and its 
chemistry is straightforward. The freeze-out process started dur- 
ing the pre-collapse phase continues up to t — 2.4 x 10 5 yr, a 
few thousand years prior to reaching point A in Fig. [3] where the 
dust temperature of 18 K results in CO evaporating again. As the 
parcel continues its inward journey and is heated up further, all 
solid CO rapidly disappears and the gas-phase abundance goes 
up to 6 x 10~ 5 at point B. During the remaining part of the infall 
trajectory, the other main carbon-bearing species are all largely 
converted into CO. At the end of the collapse (point D), 99.8% 
of all carbon and 44% of all oxygen are locked up in CO. 

The second most abundant carbon-bearing ice at the onset 
of collapse is CH4, at 1 x 10~ 5 with respect to «h- The gas- 
phase abundance of CH4 begins at 4 x 10~ 9 , about a factor of 
2500 lower. At point A, the evaporation of CO provides the first 
increase in *(CH4) through a chain of reactions starting with 
the formation of C + from CO and He + . The successive hydro- 
genation of C + produces CH^, which reacts with another CO 



molecule to form CH4: 

CO + He + -> C + + O + He, (33) 

C + + H 2 -> CH+ + hv, (34) 

CH+ + H 2 -> CH3 + H, (35) 

CHJ + H 2 -> CH+, (36) 

CH+ + CO -> CH 4 + HCO + . (37) 



The CH4 ice evaporates at point B, bringing the gas-phase abun- 
dance up to 1 x 10~ 5 . So far, the abundances of CH4 and CO are 
well coupled. The link is broken when the parcel reaches point 
C, where CH4 is photodissociated, but CO is not. This difference 
arises from the fact that CO can only be dissociated by photons 
sh ortwards of 1076 A, while CH4 can be dissociated out to 1450 
A (IVisser et alJl2009ah . The 5300 K blackbody spectrum emit- 
ted by the protostar at this time is not powerful enough at short 
wavelengths to cause significant photodissociation of CO. CH4, 
on the other hand, is quickly destroyed. Its final abundance at 
point D is 6 x 10~ 10 . 

Neutral and ionised carbon show the same trends in their 
abundance profiles, with the former always more abundant by 



a few per cent to a few orders of magnitude. Both start the col- 
lapse phase at ~10~ 8 relative to «h- The increase in total density 
at 2.1 x 10 5 yr speeds up the destruction reactions (mainly by 
OH and O2 for C and by OH and H2 for C + ), so the abundances 
go down to x(C) = 4 x 10~ 1() and x(C + ) = 5 x 10" 11 just outside 
point A. This is where CO begins to evaporate, and as a result, 
the C and C + abundances increase again. As the parcel continues 
to fall in towards point B, the evaporation of O2 and NO and the 
increasing total density cause a second drop in C and C + . Once 
again, though, the drop is of a temporary nature. Moving on to- 
wards point C, the parcel gets exposed to the stellar UV field. 
The photodissociation of CH4 leads - via intermediate CH, CH2 
or CH3 - to neutral C, part of which is ionised to also increase the 
C + abundance. Finally, at point D, the photoprocesses no longer 
play a role, so the C and C + abundances go back down. Their 
final values relative to «h are 7 x 10~ 12 and ~10~ 14 . 

5.1.3. Nitrogen chemistry 

The most common nitrogen-bearing species at t — is solid N2, 
with an abundance of 5 x 10~ 6 . The corresponding gas-phase 
abundance is 4 x 10~ 8 . The evolution of N2 parallels that of CO, 
because they hav e similar binding ene rgies and are both very 
stable molecules (Bissc hop et al.ll2006h . N2 continues to freeze 
out slowly until it gets near point A in Fig. [3] where the grain 
temperature of 18 K causes all N2 ice to evaporate. The gas- 
phase N2 remains intact along the rest of the infall trajectory 
and its final abundance is 1 x 10~ 5 , accounting for 77% of all 
nitrogen. 

The second largest nitrogen reservoir at the onset of collapse 
is NH3, with solid and gas-phase abundances of 8 x 10~ 6 and 
2 x 10~ 8 . The gas-phase abundance receives a short boost at point 
A due to the evaporation of N2, followed by 



N 2 + He + -> N + + N + He, (38) 

N + + H 2 -> NH + + H, (39) 

NH + + H 2 -> NH 2 + H, (40) 

NH 2 + H 2 -> NH3 + H, (41) 

NH3 + H 2 -> NH; + H, (42) 

NH+ + e" -> NH 3 + H . (43) 



The binding energy of NH3 is intermediate to that of O2 and 
H2O, so it evaporates between points B and C. Like H2O, NH3 
is photodissociated upon evaporation. As the last of the NH3 ice 
leaves the grains at R — 50 AU (10 AU outside of point C), the 
gas-phase reservoir is no longer replenished and x(NH3) drops to 
~10~ 14 . Some NH3 is eventually reformed as the parcel gets into 
the disk, and the final abundance at point D is 1 x 10 -10 relative 
to n H - 

With an abundance of 6 x 10~ 6 , solid NO is the third major 
initial nitrogen reservoir. Gaseous NO is a factor of twenty less 
abundant at t — 0: 3 x 10~ 7 . The NO gas is gradually destroyed 
prior to reaching point A by continued freeze-out and reactions 
with H + and H3 . It experiences a brief gain at point A from the 
evaporation of OH and its subsequent reaction with N to give 
NO and H. As the parcel continues to point B, the solid NO 
begins to evaporate and the gas-phase abundance rises to 6 x 
10 . Photodissociation reactions then set in around R = 100 
AU and the NO abundance goes back down to 6 x 10~ 9 . The 
evaporation and photodissociation of NH3 cause a brief spike in 
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the NO abundance through the reactions 

NH 3 + hv -> NH 2 + H , (44) 
NH 2 + O -> HNO + H , (45) 
HNO + hv —> NO + H. (46) 

The evaporation of the last of the NH3 ice at R — 50 AU 
eliminates this channel and the NO gas abundance decreases to 
3 x 10~ 9 at point C. NO is now mainly sustained by the reac- 
tion between OH and N. As described above, the OH abundance 
drops sharply at R = 17 AU, and the NO abundance follows suit. 
The final abundance at point D is ~1(T 14 . 

The last nitrogen-bearing species from Fig. [2] is atomic N 
itself. It starts at an abundance of 1 x 10~ 7 and slowly freezes 
out to an abundance of 2 x 10~ 8 just before reaching point A. 
At point A, N2 evaporates and is partially converted to N 2 H + . 
The dissociative recombination of N 2 H + mostly reforms N 2 , but 
it also produces some NH and N. The N abundance jumps back 
to 1 x 10~ 7 and remains nearly constant at that value until the 
parcel reaches point B, where NO evaporates and reacts with N 
to produce N 2 and O. This reduces x(N) to a low of 5 x 1CT 10 be- 
tween points B and C. Moving in further, the parcel gets exposed 
to the stellar UV field, and NO and NH3 are photodissociated to 
bring the N abundance to a final value of 5 x 10~ 6 relative to 
«h- As such, it accounts for 22% of all nitrogen at the end of the 
collapse. 

5.2. Other parcels 

At the end of the collapse (f = f acc = 2.52 x 10 5 yr), the parcel 
from Sect. 15.11 (hereafter called our reference parcel) is located 
at R = 6.3 AU and z = 2.4 AU, about 0.2 AU below the sur- 
face of the disk. As shown in Fig. [3] its trajectory passes close to 
the outflow wall, through a region of low extinction. This results 
in the photodissociation or photoionisation of many species. At 
the same time, the parcel experiences dust temperatures of up to 
150 K (Fig. O, well above the evaporation temperature of H 2 
and all other non-refractory species in our network. Material that 
ends up in other parts of the disk encounters different physical 
conditions during the collapse and therefore undergoes a differ- 
ent chemical evolution. This subsection shows how the absence 
or presence of some key chemical processes, related to certain 
physical conditions, affects the chemical history of the entire 
disk. Table [3] lists the abundances of selected species at four 
points in the disk at f acc . 

5.2.1. Oxygen chemistry 

The main oxygen reservoir at the onset of collapse is H 2 ice 
(Sect. [4]). Its abundance remains constant at 1 x 10~ 4 in our refer- 
ence parcel until it gets to point C in Fig. [3] where it evaporates 
from the dust and is immediately photodissociated. When the 
parcel enters the disk, some H 2 is reformed to produce final 
gas-phase and solid abundances of ~10" 8 relative to «h (Sect. 

ELD. 

Figure |4] shows the disk at f acc , divided into seven zones ac- 
cording to different chemical evolutionary schemes for H 2 0. The 
fraction of the disk mass in each zone is indicated. The material 
in zone 1 (61% of the disk mass) is the only material in the disk 
in which H 2 never evaporates during the collapse, because the 
temperature never gets high enough. The H 2 ice abundance is 
constant throughout zone 1 at f acc at ~1 x 10~ 4 . For the material 
ending up in the other six zones, H 2 evaporates at some point 
during the collapse phase. Zone 2 (1.0%) contains our reference 



Table 3. Abundances of selected species (relative to «h) and 
physical conditions at t — f acc at four positions in the disk (two 
on the midplane, two at the surface). 



Species 


R = 6 AU 


R = 24 AU 


R = 6 AU 


R = 24 AU 




z = 0AU 


z = AU 


z = 2 AU 


z = 6 AU 


H 2 


5(-7) 


1(-12) 


l(-8) 


4(-9) 


H 2 ice 


l(-4) 


l(-4) 


7(-8) 


8(-5) 


2 


K-6) 


6(-8) 


K-7) 


l(-5) 


O2 ice 


<1(-12) 


<1(-12) 


<1(-12) 


<1(-12) 


O 


K-6) 


l(-6) 


K-4) 


8(-7) 


OH 


2(-ll) 


<1(-12) 


<1(-12) 


9(-10) 


CO 


71-5) 


7(-5) 


8(-5) 


6(-5) 


CO ice 


<1(-12) 


<1(-12) 


<1(-12) 


<1(-12) 


CH 4 


71-8) 


l(-7) 


2C-9) 


H-5) 


CH 4 ice 


<1(-12) 


<1(-12) 


<1(-12) 


<1(-12) 


HCO+ 


1(-12) 


5(-ll) 


3(-12) 


4(-12) 


C 


<1(-12) 


<1(-12) 


9(-12) 


<1(-12) 


C + 


<1(-12) 


<1(-12) 


<1(-12) 


<1(-12) 


N 2 


l(-5) 


l(-5) 


l(-5) 


5(-6) 


N2 ice 


<1(-12) 


<1(-12) 


<1(-12) 


<1(-12) 


NH 3 


5(-10) 


2(-10) 


4(-10) 


8(-6) 


NH 3 ice 


<1(-12) 


5(-9) 


<1(-12) 


5(-8) 


NO 


l(-8) 


5(-10) 


<1(-12) 


6(-6) 


NO ice 


<1(-12) 


<1(-12) 


<1(-12) 


<1(-12) 


N 2 H + 


<1(-12) 


<1(-12) 


<1(-12) 


<1(-12) 


N 


5(-H) 


l(-9) 


5(-6) 


4(-ll) 


»h (cnr J ) 


5(12) 


8(11) 


5(11) 


4(10) 


7d(K) 


107 


76 


96 


74 


H 2 zone" 


7 


7 


2 


1 


CH 4 zone" 


5 


5 


2 


1 


NH 3 zone" 


8 


8 


5 


2 



Notes. <n) The zones from Figs. |4j [7] and |9] in which each position is 
located. 




Fig. 4. Schematic view of the history of H 2 gas and ice 
throughout the disk. The main oxygen reservoir at f acc is indi- 
cated for each zone; the histories are described in the text. The 
percentages indicate the fraction of the disk mass contained in 
each zone. Note the disproportionality of the R and z axes. The 
colours have no specific meaning other than to distinguish the 
different zones. 



parcel, so its H 2 history has already been described. The total 
H 2 abundance (gas and ice combined) at f acc is ~10~ 8 at the top 
of zone 2 and ~10~ 6 at the bottom. The gas-ice ratio goes from 
~1 at the top to ~10~ 6 at the bottom. 

The H 2 history of zone 3 (4.7%) is the same as that of 
zone 2, except that it finishes with a gas-ice ratio larger than 
unity. In both cases H 2 evaporates and is photodissociated prior 
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Fig. 5. Infall trajectory for a parcel ending up in zone 6 from Fig. 
|U showing where H2O is predominantly present as ice (black) 
or gas (grey). The two diamonds mark the time in years after the 
onset of collapse. 
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Fig. 6. Qualitative evolution of some abundances towards the 
seven zones with different H2O histories from Fig. [4] The hori- 
zontal axes show the time (increasing from left to right) and are 
non-linear. The position of points A, B, C and D from Fig. [3]is 
indicated for zone 2, which contains our reference parcel. 



to entering the disk (point C in Fig. [3}, and it partially reforms 
inside the disk (point D). Parcels ending up in zone 4 (0.4%) 
also experience the evaporation and photodissociation of H2O. 
However, the low extinction against the stellar radiation in zone 
4 prevents H2O from reforming like it does in zones 2 and 3. 

The material in zone 5 (0.9%) has a different history from 
that in zones 2—4 because it enters the disk earlier: between 1 .3 x 
10 5 and 2.3 x 10 5 yr. The material in zones 2-4 all accretes after 
2.4 X 10 5 yr. The infall trajectories terminating in zone 5 do not 
pass close enough to the outflow wall or the inner disk surface 
for photoprocesses to play a role. All H2O in zone 5 is in the gas 
phase at f acc (abundance: 1 x 10~ 4 ) because it lies inside the disk's 
snow line. The evaporation of H2O ice does not occur until the 
material actually crosses the snow line. Prior to that point, the 
temperature never gets high enough for H2O to leave the grains. 

The material ending up in zones 6 and 7 (1.6 and 30%) ac- 
cretes even earlier than that ending up in zone 5: at t — 4x 10 4 yr. 
The disk at that time is only 2 AU large and several 100 K hot, so 
the ice mantles are completely removed. Upon reaching the disk, 
the zone 6 and 7 material is swept up as the disk expands to con- 
serve angular momentum, resulting in trajectories like the one 
in Fig. [5] Because the zone 5 material enters the disk at a later 
time, its trajectories do not feature a strong outward component 
and terminate at smaller radii than the zone 6 and 7 trajectories. 
Along the outward part of the latter, material cools down and 
H 2 returns to the solid phase (Fig. |5J- At t = 1.7 x 10 5 yr, 
the parcel starts moving inwards again and comes close enough 
to the protostar for H2O to evaporate a second time. All parcels 
ending up in zone 6 have similar trajectories and the same qual- 
itative H2O history. The parcels ending up in zone 7 also have 
H2O evaporating during the initial infall and freezing out again 
during the outward part of the trajectory. However, they do not 
terminate close enough to the protostar for H2O to desorb a sec- 
ond time. Therefore, most H2O in zone 7 is on the grains at f aC c- 

Our model is not the first one in which par t of the disk 
contai n s H2O that evapo rate d and readsorbed . iLunine et alj 
(ll99ll) . IOwen et alJ (1 19921) and lOwen & Bar-Nunl (1 19931) argued 
that the accretion shock at the surface of the disk is strong 
enough for H2O to evaporate . Howeve r, based on the model of 
iNeufeld & Hollenbachl d 19941) . iPaper II showed that most of the 
disk material does not pass through a shock that heats the dust 
to 100 K or more. Moreover, the material that does get shock- 
heated to that temperature accretes close enough to the star that 
the stellar radiation already heats it to more than 100 K. Hence, 
including the accretion shock explicitly in our model would at 
most result in minor changes to the chemistry. 



As discussed in Sect. 15.1.11 H2O controls part of the oxygen 
chemistry along the infall trajectory of our reference parcel, and 
it does so for other parcels, too. Figure [6] presents a schematic 
view of the chemical evolution of six oxygen-bearing species 
towards each of the seven zones from Fig. [4] The abundances 
are indicated qualitatively as high, intermediate or low. The hor- 
izontal axes (time, increasing from left to right) are non-linear 
and only indicate the order in which various events take place. 

For material ending up in zone 1, H2O never evaporates, but 
O2 does. OH is initi tally relatively abundant but most of it disap- 
pears when the overall density increases and the reactions with 
O, N and H3 become faster. The abundance of atomic O expe- 
riences a drop at the same time, but it quickly goes back up due 
to the evaporation of CO, O2 and NO, followed by Reactions 
(l30t-(l32b. O2 also evaporates on its way to zones 2, 3 and 4, 
and because it passes through an area of low extinction, it is sub- 
sequently photodissociated. This, together with the photodisso- 
ciation of OH and H2O, causes an increase in the abundance of 
atomic O. Zones 2 and 3 are sufficiently shielded against the stel- 
lar UV field, so O2 is reformed at the end of the trajectory. This 
does not happen in the less extincted zone 4. 

Material ending up in zone 5 has the same qualitative history 
for gas-phase and solid O2 as has material ending up in zone 1 . 
The evolution of atomic O initially also shows the same pattern, 
but it experiences a second drop as the total density becomes 
even higher than it does for zone 1. The higher density also 
causes an additional drop in the OH abundance. En route towards 
zones 6 and 7, O2 evaporates and is photodissociated during the 
early accretion onto the small disk. It is reformed during the out- 
ward part of the trajectory and survives until f acc . Atomic O is 
also relatively abundant along most of the trajectory, although 
always one or two orders of magnitude below O2. In zone 6, the 
O abundance decreases at the end due to the reaction with evap- 
orating CH3, producing H2CO and H. Zone 7 does not get warm 
enough for CH3 to evaporate, so O remains intact. 

5.2.2. Carbon chemistry 

The two main carbon reservoirs at the onset of collapse are CO 
ice and CH4 ice (Sect. HJ. Their binding energies are relatively 
low (855 and 1080 K), so they evaporate throughout the core 
soon after the collapse begins (Sect. 15.1.21 ). For our reference 
parcel, the main difference between the evolution of CO and 
CH4 is the photodissociation of the latter (near point C in Fig. 
while the former remains intact. Following the analysis for 
H2O, the disk at f acc could be divided into several zones accord- 
ing to different chemical evolution scenarios for CO. However, 
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Fig. 7. As Fig.|4j but for CH4 gas and ice. The main carbon reser- Fig. 9. As Fig. |4j but for NH3 gas and ice. The main nitrogen 
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Fig. 8. As Fig. [6] but for the five zones with different CH4 histo- 
ries from Fig. [7] 



the entire disk has the same qualitative CO history: CO does 
not undergo any processing apart from evaporating early on in 
the collaps e phase, and it is the main carbon reservoir through- 
out (>99%; Pap er]}) ■ Hence, the disk is divided according to the 
evolution of CH4 instead (Fig. [7j». 

For material that ends up in zone 1 (representing 62% of the 
disk mass), the only major chemical event for CH4 is the evapo- 
ration during the initial warm-up of the core. It is not photodis- 
sociated at any point during the collapse, nor does it freeze out 
again or react significantly with other species. Material ending 
up closer to the star, in zones 2 and 3 (5.3 and 0.4%), is suffi- 
ciently irradiated by the stellar UV field for CH4 to be photodis- 
sociated. The amount of extinction increases towards the end of 
the trajectories terminating in these zones, allowing some CH4 
to reform. The final abundance ranges from 10~ 10 in zone 3 to 
10~ 7 in the most shielded parts of zone 2. 

Zone 4 (0.9%) contains material that accretes onto the disk 
several 10 4 yr earlier than does the material in zones 2 and 3. It 
always remains well shielded from the stellar radiation, so the 
only processing of CH4 is the evaporation during the early parts 
of the collapse. The CH4 history of zones 1 and 4 is thus qualita- 
tively the same. Zone 5 (32%) consists of material that accretes 
around t — 4x 10 4 yr and is subsequently transported outwards to 
conserve angular momentum (Fig. [5). CH4 in this material evap- 
orates before entering the young disk and is photodissociated as 
it gets within a few AU of the protostar. The resulting atomic C 
is mostly converted into CO and remains in that form for the rest 
of the trajectory. Hence, even though the extinction decreases 
again when the parcel moves outwards, no CH4 is reformed. 

The evolution of the abundances of CH4 gas and ice, CO gas 
and ice, HCO + , C and C + towards each of the five zones is plot- 



ted schematically in Fig. [8] The HCO + abundance follows the 
CO abundance at a ratio that is roughly inversely proportional to 
the overall density. Hence, the HCO + evolution is qualitatively 
the same towards each zone: it reaches a maximum abundance 
of a few 10~ 10 when CO evaporates and gradually disappears 
as the density increases along the rest of the infall trajectories. 
The most complex history amongst these seven carbon-bearing 
species is found in C and C + . Towards all five zones, they are ini- 
tially destroyed by reactions with H2, O2 and OH. Some C and 
C + is reformed when CO evaporates, but the subsequent evapo- 
ration of O2 and OH causes the abundances to decrease again. 
En route to zones 2, 3 and 5, the photodissociation of CH4 leads 
to a second increase in C and C + , followed by a third and final 
decrease when the parcel moves into a more shielded area. 

5.2.3. Nitrogen chemistry 

Most nitrogen at the onset of collapse is in the form of solid 
N 2 (41%), solid NH 3 (32%) and solid NO (22%). The evolution 
of N2 during the collapse is the same as that of CO, except for 
a minor difference in the binding energy. Both species evapo- 
rate shortly after the collapse begins and remain in the gas phase 
throughout the rest of the simulation. Neither one is photodisso- 
ciated, because the protostar is not hot enough to provide suffi- 
cient UV flux shortwards of 1 100 A. 

The evolution of NH3 shows a lot more variation, as illus- 
trated in Fig. [9] The disk at f acc is divided into eight zones with 
different NH3 histories. No processing occurs towards zone 1 
(8.3% of the disk mass): the temperature never exceeds the 73 
K required for NH3 to evaporate, so it simply remains on the 
grains the whole time. Material ending up in zone 2 (25%) does 
get heated above 73 K, and NH3 evaporates. However, it freezes 
out again at the end of the trajectory because zone 2 itself is not 
warm enough to sustain gaseous NH3. The final solid NH3 abun- 
dance in zones 1 and 2 is about 8 x 10~ 5 relative to «h (TableQ. 

NH3 ending up in zone 3 (20%) evaporates from the grains 
just before entering the disk, but it is immediately photodissoci- 
ated or destroyed by HCO + . Towards the end of the trajectory, 
some NH3 is reformed from the dissociative recombination of 
NH^. It rapidly freezes out to produce a final solid NH3 abun- 
dance of 10 -10 — 10 8 . Zone 4 (7.2%) has the same history, except 
that there is an additional adsorption-desorption cycle before the 
destruction by photons and HCO + . 

Our standard parcel ends up in zone 5 (7.7%); its NH3 evolu- 
tion is mainly characterised by photodissociation above the disk 
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Fig. 10. As Fig. [6] but for the eight zones with different NH3 
histories from Fig. [9] 



and reformation inside it (Sect. 15.1.3b . NH3 evaporates when it 
gets to within about 200 AU of the star, halfway between points 
B and C in Fig. [3] It is immediately photodissociated, but some 
NH3 is reformed once the material is shielded from the stellar 
radiation. The reformed NH3 remains in the gas phase. No ref- 
ormation takes place in the less extincted zone 6 (0.4%), which 
otherwise has the same NH3 history as zone 5. 

Material ending up in zone 7 (0.9%) does not pass close 
enough to the outflow wall for NH3 to be photodissociated upon 
evaporation. Instead, NH3 is mainly destroyed by HCO + and at- 
tains a final abundance of ~10~ 9 . Lastly, zone 8 (30%) contains 
again the material that accretes onto the disk at an early time and 
then moves outwards to conserve angular momentum. Its NH3 
evaporates already before reaching the disk and is subequently 
dissociated by the stellar UV field. As the material moves away 
from the star and is shielded from its radiation, some NH3 is 
reformed out of NH4 to a final abundance of ~10~ 10 . 

The abundances of N2FF and atomic N are largely controlled 
by the evolution of N2 and NH3, as shown schematically in 
Fig. [TO] In all parcels, regardless of where they end up, N2I-F 
is mainly formed out of N2 and Hi, so its abundance goes up 
when N2 evaporates shortly after the onset of collapse. It gradu- 
ally disappears again as the collapse proceeds due to the increas- 
ing density. The atomic N abundance at the onset of collapse is 
1 x 10~ 7 . It increases when N2 evaporates and decreases again a 
short while later when NO evaporates (Sect. lBTOI l. For material 
that ends up in zones 1 and 2, the N abundance is mostly con- 
stant for the rest of the collapse at a value of 10~ 13 (inner part of 
zone 2) to 10~ 10 (outer part of zone 1). Material that ends up in 
the other six zones is exposed to enough UV radiation for NO 
and NH3 to be photodissociated, so there is a second increase in 
atomic N. Zones 3, 4, 7 and 8 are sufficiently shielded at f acc to 
reform some or all NO and NH3, and the N abundance finishes 
low. Less reformation is possible in zones 5 and 6, so they have 
a relatively large amount of atomic N at the end of the collapse 
phase. 



5.2.4. Mixing 

Given the dynamic nature of circumstellar disks, the zonal dis- 
tribution presented in the preceding subsections may offer too 
simple a picture of the chemical composition. For example, there 
are as yet no first principles calculations of the processes respon- 
sible for the viscous tran sport in disks. The radial ve locity equa- 
tion used in our model (I Visser & Dullemondl l2010) is suitable 
as a zeroth-order description, but it omits the effects of m ixing 
dLewis & Prind [19801 : ISemenov et al. 2006; lAikawall2007l ) and 
it cannot explain important observational features like episodic 



accret ion dKenvon & Hartmannll 19951) . Indeed, with the lToomrel 
( 1964) Q parameter reaching a minimum of 1 .4 at f acc , the disk is 
only barely gravitationally stable. Small changes in the model or 
in the initial conditions would suffice to make the disk unstable 
at various points during the collapse phase, resulting in a redistri- 
bution of material beyond the equations currently used. Hence, 
in real disks, both the shapes and the locations of the zones are 
likely to differ from what is shown in Figs. [HQ and [9] A larger 
degree of mixing would also make the borders between the zones 
more diffuse than they are in our simple schematic representa- 
tion. Nevertheless, the general picture from this section offers a 
plausible description of the chemical history towards different 
parts of the disk. Spectroscopic observations at 5-10 AU resolu- 
tion, as possible with ALMA, are required to determine to what 
extent this picture holds in reality. Based on the abundances in 
Table[3] potential tracers are H 2 0, NO and NH 3 . 



6. Chemical history versus local chemistry 

Section [5] contains many examples of abundances increasing or 
decreasing on short timescales of less than a hundred years (see, 
e.g., Fig. [2|. It appears that the abundances respond rapidly to 
the changing physical conditions as material falls in superson- 
ically through the inner envelope and accretes onto the disk. 
However, this does not necessarily mean that the abundances are 
always in equilibrium. The current section explores the ques- 
tion of whether the disk is in chemical equilibrium at the end 
of the collapse, or whether its chemical composition is a non- 
equilibrium solution to the conditions encountered during the 
collapse. To that end, the chemistry is evolved for an additional 
1 Myr beyond f acc . The density, temperature, UV flux and ex- 
tinction are kept constant at the values they have at ? acc , and all 
parcels of material are kept at the same position. Clearly, this is 
a purely hypothetical scenario. In reality, the disk would change 
in many ways after f acc : it spreads in size, the surface layers be- 
come more strongly irradiated, the temperature changes, the dust 
coagulates into planetesimals, gas is photoevaporated from the 
surface layers, and so on. All of these processes have the po- 
tential to affect the chemical composition. However, they would 
also interfere with the attempt to determine whether the disk is 
in chemical equilibrium at f acc . This question is most easily ad- 
dressed by evolving the chemistry for an additional period of 
time at constant conditions. Hence, all of the physical changes 
that would occur during the post-collapse phase are ignored. 

The 21 species discussed in Sect. [5]can be divided into two 
categories: those whose abundance profile changes during the 
post-collapse phase, and those whose abundance profile remains 
practically the same. Members of the "changed" category are 
O2, OH, CH 4 , NH 3 and NO (all gaseous), NH 3 ice, and atomic 
O and N. The thirteen species in the "unchanged" category are 
H2O gas and ice, O2 ice, CO gas and ice, CH4 ice, C, C + , N2 
gas and ice, NO ice, HCO + and N2H + . The individual gas and 
ice abundances are summed in Fig. [TT] and compared at t — f acc 
and f acc + 1 Myr. The total H2O, CO and N2 abundances do not 
change significantly during the post-collapse phase, while the 
total O2, CH4, NH3 and NO abundances change by more than 
two orders of magnitude in a large part of the disk. 

There are two areas in the disk where the abundances of all 
21 species remain nearly constant: near the surface out to ~10 
AU, and at the midplane between ~5 and ~25 AU. The chem- 
istry in the first area is dominated by fast photoprocesses. The 
second area, near the midplane, is the densest part of the disk 
and therefore has high collision frequencies. In both cases, the 
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Fig. 11. Abundances of total H 2 0, 2 , CO, CH 4 , N 2 , NH 3 and 
NO (gas and ice combined) throughout the disk at the end of 
the collapse phase (f = f acc ) and after an additional 1 Myr post- 
collapse phase (f = f acc + 1 Myr). 

chemical timescales are short and the chemistry reaches equilib- 
rium during the final stages of accretion. 

Outside these two "equilibrium areas", one of the key pro- 
cesses during the post-collapse phase is the conversion of gas- 
phase oxygen-bearing species into H 2 0, which subsequently 
freezes out. Because H 2 ice is already one of the main oxy- 
gen reservoirs at f acc throughout most of the disk (Fig. @J, its 
abundance only increases by about 20% during the rest of the 
simulation. These 20% represent 10% of the total oxygen budget 
and correspond to 2 x 10~ 5 oxygen nuclei per hydrogen nucleus. 

One of the species converted into H 2 ice is NO, which re- 
acts with NH 2 through 



NH 2 + NO -» N 2 + H 2 0, 



(47) 



followed by adsorption of H 2 onto dust. This reaction is re- 
sponsible for the post-collapse destruction of gas-phase NO out- 
side the two "equilibrium areas", as well as for a ~20% increase 
in gas-phase N 2 . NO is the main destructor of atomic N via 



NO + N -> N 2 + O, 



(48) 



so the conversion of NO into H 2 leads to a higher N abundance. 
The NH 2 required for Reaction ( |47| ) is formed from NH3: 



NH, + HCO + 



NH^ + CO, 



(49) 

NH4 + e" -> NH 2 + H 2 /2H. (50) 

The gas-phase reservoir of NH3 is continously fed by evapo- 
ration of NH 3 ice, so Reactions (|49l ), (T50b and ( |4"7T i effectively 
transform all NH3 into N 2 . 

Two other important species that are destroyed in the post- 
collapse phase are gaseous 2 and CH4. 2 is ionised by cosmic 
rays to produce 2 , which reacts with CH4: 



CH 4 + OJ 
HCOOH+ + e 



HCOOH^ + H, 
HCOOH + H . 



(51) 
(52) 



The formic acid (HCOOH) freezes out, thus acting as a sink for 
both carbon and oxygen. At t — f acc + 1 Myr, the solid HCOOH 
abundance is 3 x 10~ 6 , accounting for 3% of all oxygen and 4% 
of all carbon. The presence of s uch a sink is a common feature 
of disk chemistry models (e.g.. lAikawa & Herbstl [l999). Gas- 
phase species are processed by He + and HJ until they form a 
species whose evaporation temperature is higher than the dust 
temperature. This is HCOOH in our case, but it could also be 
carbon-chain molecules like C 2 H 2 or C3H4. 

Another important post-collapse sink reaction in our model 
is the freeze-out of HNO. At t > f acc , OH is primarily formed by 



HNO + O 



NO + OH 



(53) 



throughout most of the disk. This reaction is also one of the main 
destruction channels for atomic O. Hence, the gradual freeze-out 
of HNO after f acc leads to a decrease in OH and an increase in O. 
Also contributing to the higher O abundance is the fact that OH 
itself is an important destructor of O through Reaction (fT3l . 

The post-collapse chemistry described here is merely meant 
as an illustration of what might take place in a real circumstellar 
disk after the main accretion phase comes to an end. Depending 
on how the disk continues to evolve physically, other reac- 
tions may become more important than the ones listed above. 
However, with the abundances of several key species changing 
by orders of magnitude in our post-collapse experiment, the disk 
was clearly not in chemical equilibrium at the end of the col- 
lapse. Even the two "equilibrium areas" (near the surface out to 
10 AU and from 5 to 25 AU at the midplane) are not truly in 
chemical equilibrium when considering less abundant species. 
For example, at 20 AU on the midplane, post-collapse process- 
ing increases the abundance of solid HCOOH by a factor of 13 
(from 7 x 10" 8 to 9x 10~ 7 relative to ra H ) and that of solid CH3OH 
by a factor of 27 (from 3 x 10~ 13 to 8 x 10~ 12 ). 

The lack of chemical equilibrium at the end of the collapse 
phase means that one has to be careful in choosing the initial 
abundances when running T Tauri or Herbig Ae/Be disk chem- 
istry models. The abundances of minor species, such as most or- 
ganic molecules, can be particularly sensitive to what happened 
before the disk was formed. Ideally, any disk chemistry model 
should therefore include a realistic core collapse phase to set the 
stage for the chemical evolution in the disk itself. 



7. Comets: mixed origins? 

The chemical composition of cometary ices shows many simi- 
larities to that of young protostars, but the abundance of a given 
species may vary by more than an order of magnitude from 
one c omet to the next (lA'Hearn et al.ll 99*51: iMumma & Charnlevl 
l201ll) . Both points are illustrated in the top panel of Fig. 
[T2l where the abundances of several species in the comets 
lP/Halley, C/1995 Ol (Hale-Bopp), C/1996 B2 (Hyakutake), 
C/1999 HI (Lee), C/1999 S4 (LINEAR) and 153P/Ikeya-Zhang 
are plotted against the gas-phase abundances in the embedded 
protostar IRAS 16293-2422 (warm inner envelope), the bipolar 
outflow L 1157, and the four hot cores W3(H 2 0), G34.3+0.15, 
Orion HC and Orion CR. Al l abundances are gathered from 
iBockelee-Morvan et all d200d |2004l) and ISchoier et all (|2002). 



The cometary data are converted from abundances relative to 
H 2 to abundances relative to «h by assuming «(H 2 0)/«h = 
5 x 10~ 5 . 

Each point in the top panel of Fig. [12] is characterised by 
the mean value from the available sources (the diamonds) and 
the total spread in measurements (the error bars). Uncertainties 
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Fig. 12. Top: abundances observed in comets versus gas-phase 
abundances observed in young protostars. Middle and bottom: 
abundances observed in comets versus total abundances (gas 
and ice) from our model at t — f acc . All abundances are rela- 
tive to «h- For the observational data points, the error bars in- 
dicate the spread between sources. For the model data, they in- 
dicate the range of values across the entire disk {middle) or the 
comet-forming region (bottom). Dotted error bars indicate the 
full range extends to less than 1CT 12 . The diagonal dotted lines 
represent hypothetical one-to-one relationships between abun- 
dances on the x and y axes. 



from individual measurements are not shown. The dotted line 
marks the one-to-one relationship where the cometary abun- 
dances equal the protostellar abundances. The data generally fol- 
low this line, suggesting that the material ending up in comets 
underwent little chemical processing during the evolution from 
protostar to main-sequence star. However, the plot also shows 
variations of at least an order of magnitude in the cometary 
abundances for CO, H 2 CO, CH 3 OH, HNC and H 2 S, as well as 
smaller variations for other species. A possible explanation for 
the different chemical compositions is that the comets formed in 
different parts of the solar nebula. If that is indeed the case, there 
must have been some degree of chemical processing between 
the protostellar stage and the time where each comet formed. 
Another point to note is that the elemental nitrogen abundance 
in comets is at least a factor of three lower than it is in protostars 
dWvckoff et al.lfl99ll) . Although the reason for this deficiency 
remains unclear, it also points at a certain degree of chemical 
processing. 

Comets are thought to have formed at the gravitational mid- 
plane of the circumsolar disk, between 5 and 30 AU from the 
young Sun dBockelee-Morvan et alj 120041) . In our model disk, 
the 5-30 AU range is almost exactly the area containing mate- 
rial that accreted at an early time (4 x 10 4 yr after the onset of 
collapse) and stayed in the disk for the remainder of the collapse 
phase (Fig. [5]). As discussed in Sect. 15.21 this material undergoes 
a larger degree of chemical processing than any other material 
in the disk, due to its accreting close to the star. However, such a 
large degree of processing is hard to reconcile with the observed 
similarities between cometary and protostellar abundances, and 
it is therefore unlikely that comets form entirely out of mate- 
rial that accreted close to the star. Instead, as argued below, our 
model results suggest that comets consist of a mix of processed 
and unprocessed material. 

The middle and botto m panels of Fig. [12] compare the 
cometary abundances from lBockelee-Morvan et alJd2000l [2004) 
to the volume-averaged abundances from our model at the end 
of the collapse phase. For all species, the gas and ice abundances 
from the model are summed and displayed as one. The horizon- 
tal error bars show again the spread in abundances between in- 
dividual comets. In the middle panel, the vertical bars show the 
spread across the entire model disk; in the bottom panel, they 
show the spread in the 5-30 AU region at the midplane. Dotted 
extensions to the error bars indicate that the range of values ex- 
tends to less than 10~ 12 . 

The disk-averaged abundances tend to cluster around the hy- 
pothetical one-to-one relationship indicated by the diagonal dot- 
ted line in the middle panel of Fig. [12] Two notable exceptions 
are CH3OH and HCOOCH3. Both species are predominantly 
formed through grain-surface chemistry, which is not included 
in our model. The spread in model abundances is at least six 
orders of magnitude, but given the wide range of physical condi- 
tions that are sampled, such a large spread is to be expected. The 
5-30 AU abundances from the bottom panel of Fig. Q~2]are lower 
for many species than the disk-wide abundances. For example, 
the average NH3 abundance decreases by more than three orders 
of magnitude. NH3 is photodissociated along the trajectories that 
terminate in the comet-forming zone, whereas other trajectories 
are more shielded (Sect. lBTOT l. 

The abundances in the bottom panel of Fig. [12] are sampled 
over a smaller region, so they show less variation than those in 
the middle panel. However, they still take on a wider range of 
values than do the cometary abundances. More importantly, the 
5-30 AU data show no correlation with the comet data. There is 
no indication from Sect. [6] that post-collapse processing brings 
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the abundances back to near what they were in the protostellar 
stage. Hence, the midplane material between 5 and 30 AU in 
our model does not appear to be analogous to the material from 
which the solar-system comets formed. 

How plausible is the large degree of chemical processing for 
material that ends up in the comet-forming zone? The amount 
of processing depends on the range of physical conditions en- 
countered along a given infall trajectory, which in turn depends 
mostly on how close to the protostar the trajectory gets. The 
physics in our model are kept as simple as possible, so a trajec- 
tory similar to the one drawn in Fig.[3]may not be fully realistic. 
On the other hand, the back-and-forth motion results from the 
well-known concept of conservation of angular momentum. As 
the inner parts of the disk accrete onto the star, the outer parts 
must move out to maintain the overall angular momentum. They 
may be pushed inwards again at a later time if a sufficient amount 
of mass is accreted from the envelope at larger radii. This hap- 
pens in our simple model, but also in the more physica lly realis- 
tic hydrodynamical simulations of lBrinch" etalJ (120081) . 

Accepting some degree of back-and-forth motion as a com- 
mon part in disk evolution, the question remains how close to 
the star the material gets before moving out to colder parts of 
the disk. Our revised solution to the problem of sub-Keplerian 
accretion dVisser & Dullemondll20lol) res ults in r adial velocity 
profiles that are different from the ones in iPaper 1 With the old 
profiles, all material that accretes inside of the snow line al- 
ways remains inside of the snow line. The comet-forming zone 
therefore did not contain any material in which H2O had evap- 
orated and readsorbed onto the grains. Such material does exist 
in the comet-forming zone with the new velocity profiles (Fig. 
[5]). However, the result depends strongly on the initial physical 
conditio ns. If the new model is run for the parameter grid from 
Paper I, there are three cases (out of eight) where material ac- 
cretes inside of the snow line and then moves out beyond it; the 
standard disk discussed throughout this work is one of them. The 
critical question is whether the point at which the radial veloc- 
ity in the disk changes sign (from inwards to outwards) lies in- 
side or outside of the snowline. In general, the radial velocity 
changes sign at smaller radii for less massive disks. Such sys- 
tems are therefore more likely to have material accreting inside 
of the snowline and then moving out. Indeed, the three cases in 
our parameter grid where this happens are the three with the least 
massive disks. Quantitatively, the dividing factor appears to be a 
ratio of 0.2 between the disk mass at t — f acc and the core mass at 
t — 0. For larger ratios, material that accretes inside of the snow 
line always remains there. 

The presence of crystalline silicate dust in disks provides 
a strong argument in favour of trajectories similar to the one 
from Fig. [5] Observations show crystalline fractions at R ss 10 
AU that are significant ly larger than what is found in the ISM 
dBouwman et al.l l2008h . Crystallisation of amorphous silicates 
requires temperatures of 800 K. The disk is less than 100 K at 
10 AU, so a possible explanation for the high crystalline frac- 
tions is that material accreted in the much hotter inner disk and 
then moved outwards. Although the particular trajectory from 
Fig. |5] does not reach the required temperature of 800 K, other 
trajectories do. Crystalline dust has also been detected in several 
comets, including C/1995 01 (Hale-Bopp) and 8 IP/Wild 2, pro- 
viding direct evidence that part of their constituent material has 
been heated t o temperatures well above the evaporati on temper- 
ature of H 2 (IWooden et al.lll999HKeller et al.ll2006l) . 

If Hale-Bopp, Wild 2 and other comets contain crystalline 
silicates, they must also contain ices that underwent a large de- 
gree of chemical processing. Likewise, the presence of amor- 



phous silicates in comets is indicative of chemically unprocessed 
ices. Hence, the material from which comets are formed must 
be of mixed origins: some of it accreted close to the star and 
was heated to high temperatures, while another part accreted at 
larger radii and remained cold. Within the context of our model, 
this requires that material from beyond R = 30 AU is radi- 
ally mixed into the comet-forming zone between 5 and 30 AU. 
Alternatively, vertical mixing may bring relatively pristine mate- 
rial from higher up in the disk down into the comet-forming zone 
at the midplane. As for the chemical variations between indi- 
vidual comets, our model shows many examples of abundances 
changing with position or with time. The variations can thus be 
explained by having the comets form at different positions in the 
circumsolar disk, or at different times during the disk's lifetime. 

8. Caveats 

One of the model's three important caveats is that it sets the 
gas temperature equal to the dust temperature, rather than cal- 
culating it explicitly. This approach is valid for the interior of 
the disk (optically thick to UV radiation), but it breaks down 
in the surface layers and in the inner parts of the envelope 
(Kamp & Dullemondl 120041) . The chemistry in these optically 
thin regions is mostly controlled by photoprocesses and ion- 
molecule reactions, neither of which depend strongly on tem- 
perature. One aspect that would be affected is the gas-phase pro- 
duction of H2O through the reactions 

O + H 2 -> OH + H, (54) 
OH + H 2 -> H 2 + H, (55) 

which have activation barriers of 3100 and 2100 K , respec- 
tively dNataraian & Rothl 119871: lAtkinson etall 120041) . Hence, 
our model probably underestimates the amount of H 2 in the 
surface of the disk. 

A higher gas temperature would also change some of the 
physics in the model, such as increasing the pressure of the disk. 
This would push up the disk-envelope boundary and possibly 
change the spatial distribution of where material accretes onto 
the disk. However, the bulk of the accretion currently takes place 
in shielded regions, so only a small fraction of the infalling ma- 
terial is expected to be affected. 

Another caveat is the adopted shape of the stellar spectrum, 
which may have larger chemical consequences than the gas tem- 
perature. Presently, the stellar radiation field is a blackbody spec- 
trum at the effect ive stell ar temperature. Our star never gets hot- 
ter than 5800 K dPaper J) , which is not enough to produce UV 
photons of sufficient energy t o dissociate CO and H 2 . M any T 
Tauri stars have a UV excess dHerbig & Goodrichll 19861) . which 
would allow CO and H 2 to be photodissociated. The higher UV 
flux would also enhance the photoionisation of atomic C, proba- 
bly resulting in C + being the dominant form of carbon between 
points C and D in Fig. [3] In turn, this would boost the produc- 
tion of carbon-chain species like C 2 S and HC3N. As soon as the 
material enters the disk and is shielded again from the stellar 
radiation, C + is converted back into C and CO. However, some 
signature of the temporary high C + abundance may survive. 

The third important caveat has to do with the dust acting as 
a possible sink for the gas-phase chemistry. In our model, once 
an atom or molecule freezes out, there are three possibilities: 
(1) it may stay on the grain if the temperature and incident UV 
flux remain low; (2) it may desorb again if the temperature or 
UV flux gets high enough; or (3) it may be hydrogenated in a 
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few select cases (Sect. lTOl ). In reality, there are at least three ad- 
ditional possibilities: (4) volatile species like CO, CH4 and N2 
can be trapped in the H2O ice matrix; (5) the atom or molecule 
in question may undergo grain-surface reactions not included in 
our model; or (6) the dust grain may grow and settle to the mid- 
plane, effectively removing the molecules in its icy mantle from 
the active chemistry elsewhere. The latter scenario has been in- 
voked to explain the low H2O abun dance in the circumstellar 
disk of DM Tau dBergin et alj|2010t) . Processes 4, 5 and 6 can 
all act as a sink of gas-phase CO and other abundant species, 
potentially diminishing their dominant role in the gas phase. A 
quantitative treatment of each process is required to assess how 
much they would alter the chemical evolution presented here. 

9. Conclusions 

This paper describes the two-dimensional chemical evolution 
during the collapse of a molecular cloud core to form a low- 
mass prot ostar and a c i rcumst ellar disk. The model, desc r ibed i n 
detail bv lVisser et all d2009bl) and lVisser & Dullemondl (l2010t) . 
consists of semi-analytical density and velocity profiles through- 
out the core and the disk. The dust temperature and the UV radi- 
ation field - both important for the chemical evolution - are cal- 
culated with a full radiative transfer method. The chemical part 
of the model features a full gas-phase network, including ad- 
sorption and desorption from dust grains, as well as basic hydro- 
genation reactions on the grain surfaces. Starting from realistic 
initial conditions, the chemical evolution is computed for parcels 
of material terminating at selected positions in the disk. Special 
attention is paid to parcels ending up in the comet-forming zone. 
Our conclusions are as follows: 

• The chemistry during the collapse phase is controlled by a 
small number of key chemical processes, each of which is 
activated by changes in the physical conditions. The evap- 
oration of CO, CH4 and H2O at approximately 18, 22 and 
100 K is one set of such key processes. Another set is the 
photodissociation of CH4 and H2O (Sect. l5Tl ). 

• At the end of the collapse phase, the disk can be divided 
into several zones with different chemical histories. The dif- 
ferent histories are related to the presence or absence of the 
aforementioned key processes along various infall trajecto- 
ries. Spectroscopic observations at high spatial resolution are 
required to determine whether this zonal division really ex- 
ists, or if it is smoothed out by mixing (Sect. 15721 ). 

• Part of the material that accretes onto the disk at early times 
is transported outwards to conserve angular momentum, and 
may remain in the disk for the rest of the collapse phase. It is 
heated to well above 100 K as it accretes close to the star, so 
H2O and all other non-refractory species evaporate from the 
grains. They freeze out again when the material cools down 
during the subsequent outward transport (Sect. [5T2J- 

• When the chemistry is evolved for an additional 1 Myr at 
fixed physical conditions after the end of the collapse phase, 
the abundances of most species change throughout the disk. 
Hence, the disk is not in chemical equilibrium at the end 
of the collapse, and one should be careful in choosing the 
initial abundances when running stand-alone disk chemistry 
models (Sect. [6]). 

• Material that ends up in the comet-forming zone undergoes a 
large degree of chemical processing, including the evapora- 
tion and readsorption of H2O and species trapped in the H2O 
ice. This is consistent with the presence of crystalline sili- 
cates in comets. However, it is inconsistent with the chem- 



ical similarities observed between comets and young proto- 
stars, which are indicative of little processing. Hence, it ap- 
pears that the cometary material is of mixed origins: part of 
it was strongly processed, and part of it was not. The chem- 
ical variations observed between individual comets suggest 
they were formed at different positions or times in the solar 
nebula. Fully pristine ices only appear in the upper and outer 
parts of this particular disk model, although turbulent mixing 
may move them into the comet-forming zone (Sect. |7]i- 
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